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Abstract 

\Q . Photons propagating in strong magnetic fields are subject to a phenomenon called the "vacuum 

birefringence" where refractive indices of two physical modes both deviate from unity and are 
1 different from each other. We compute the vacuum polarization tensor of a photon in a static 

t— i ■ and homogeneous magnetic field by utilizing Schwinger's proper-time method, and obtain a series 

£> ■ 

. !_< . representation as a result of double integrals analytically performed with respect to proper-time 

X- 

^ ' variables. The outcome is expressed in terms of an infinite sum of known functions which is 

plausibly interpreted as summation over all the Landau levels of fermions. Each contribution from 
infinitely many Landau levels yields a kinematical condition above which the contribution has an 
imaginary part. This indicates decay of a sufficiently energetic photon into a fermion-antifermion 
pair with corresponding Landau level indices. Since we do not resort to any approximation, our 
result is applicable to the calculation of refractive indices in the whole kinematical region of a 
photon momentum and in any magnitude of the external magnetic field. 
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I. INTRODUCTION 



It has been long studied that structure of the quantum vacuum in QED would be modified 
in the presence of externally applied strong electromagnetic fields |2j, and modification of 
the vacuum could entail novel phenomena such as vacuum birefringence of a photon, photon 
decay into an electron-positron pair, and photon splitting. It is natural to expect these effects 
to occur because the vacuum in QED is filled with electrons in the Dirac sea, and they react 
as 'media' like in ordinary substances in response to external fields Whereas any such 
effect has not been established in experiments, recent years have witnessed an increasing 
interest in possibilities that extremely strong electromagnetic fields would be realized in 
several different situations (see, for example, Ref. {4] for a wide range of physics relatec 
to strong fields). Primary examples are relativistic (non-central) heavy- ion collisions |s — 8 1 
and strong.y-magneti.ed compact stars such as magnetars Q, both of which are thought to 
accompany electromagnetic fields far stronger than so-called the "critical field" for electronsjj 
B c = E c = m 2 J\e\. Many observable effects are proposed in relation to magnetars [u| and 
also to heavy-ion collisions s), 12]. Besides, intensity of high-field laser has been rapidly 



increasing and is now about to reach the critical field strength. Many authors have been 
addressing theoretical aspects of the possible vacuum physics that could be studied by the 
ultra-high-intensity laser (see for example, Ref. 13[). Thus, in near future, we would be able 
to study the vacuum physics in strong fields in somewhat controllable environments. With 
these current situations in mind, we should further pursue towards deeper understanding of 
the vacuum physics in strong fields. 

In the present paper, we provide a theoretical framework necessary for the description 
of vacuum birefringence and photon decay in strong magnetic fields which are typical and 
important examples of the vacuum physics in strong fields. We give an analytic expression 
for the vacuum polarization tensor of a photon within one loop of a 'dressed' fermion that 
includes all-order interaction with the external field at tree level. Notice that the strong fields 
compensate the smallness of the coupling constant, making a particular kind of higher-order 
diagrams enhanced. Moreover, when the fields are stronger than the critical ones, naive 
perturbative expansion breaks down and thus one has to sum up all-order diagrams that are 
1 Since the magnetic axis of a magnetar is generally not the same as the rotation axis (thus, a pulsar), 

the magnetic field shows strong time dependence and will induce an electric field of the similar order of 

strength. 



enhanced by the external fields. The dressed fermion propagator which appears in the one- 
loop diagram is obtained after such resummation with respect to the external fields. As a 
result of resummation of higher-order diagrams, observables acquire nonlinear dependences 
on the external field. Since interaction with the external field is encoded as a linear coupling 
at the Lagrangian level, strong-field physics describing this kind of nonlinear effect is called 
the "nonlinear QED" . 

In spite of a simple structure of the one-loop diagram, complete description of the vacuum 
polarization tensor in the whole kinematical region has not been available so far. While the 
vacuum polarization tensor at the one-loop 



respect to two proper-time variables 
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evel is g iven in an integral representation with 
L7, 3-21], its integrand is composed of rapidly 



oscillating exponential factors which have prevented previous studies from precise analytic 
and even numerical understanding of the phenomena. Because of this difficulty, our under- 
standing of the vacuum birefringence was limited to a restricted phase space of the photon 
momentum 16l4l8j] or to a strong-filed case where lowest Landau level approximation is 
applicable [ji]]. We show however that the double proper-time integral can be analytically 
performed owing to a double infinite series expansions of the integrand. Since our calcula- 
tion does not resort to any approximation, the result is applicable to any value of a photon 
momentum and any magnitude of the external field. Namely, our result can also treat, as 
well as the vacuum birefringence, the photon decay into a fermion- ant if ermion pair which 
occurs when the photon energy is sufficiently large. 

This paper is organized as follows. First, in Sec. Ull we present the integral representa- 
tion of the vacuum polarization tensor of a photon in the external magnetic field. Then, in 
Sec. IHIl we explain how the vacuum birefringence appears in strong magnetic fields. An- 
alytic evaluation of the integral is discussed in detail in Sec. IIVI In this section, we also 
provide physical interpretation of the double infinite series expansion, and then carefully 
investigate singularity structure of the vacuum polarization tensor. Kinematical condition 
of the photon decay also appears here. Summary and prospects are given in the last section 
where we mention another important step to obtain the refractive indices. In Appendices, we 
explain supplemental materials such as details of the proper-time method, renormalization 
issues, and so on. In principle, we are able to calculate the refractive indices by using the 
polarization tensor obtained in the present paper. Taking care of the procedure mentioned 



in the last section, we will show the results of refractive indices in the next paper 



22] 



II. VACUUM POLARIZATION TENSOR IN EXTERNAL FIELDS 



In this section, we provide theoretical framework for the vacuum polarization tensor of 
a propagating photon in a strong external field, which is necessary for the calculation of 
the refractive indices. In particular, we discuss the case of a strong magnetic field, but the 
techniques developed in this section should be equally applied to a case with a strong electric 
field, which will be reported separately. 

We consider a system of massive charged fermions interacting with photons. Thus we 
start with the standard spinor QED Lagrangian: 

JSf = i> (0 - m) t/; - -F^F^ , (1) 



where we adopt a convention for the covariant derivative, D M = <9 M + ieA^ix) |23|]. Here, 
"e" and "m" representatively denote charge and mass of a fermion (e is negative for an 
electron). Since the gauge field contains both the external field and dynamical 
(i.e., propagating) degrees of freedom a M : A^ = A^ + a M , the fermion kinetic term in 
the Lagrangian ([I]) determines the coupling of the fermion to the external field. Whereas 
the fundamental QED Lagrangian (CQ) describes only the linear interaction at the classical 
level, some of higher-order quantum effects become important when they are enhanced 
due to the strong external field, and give rise to nonlinear interaction among photons. 
For example, when the external field is a strong magnetic field, insertion of an external 
field line gives a factor of 0(\eB\/m 2 ) and n external field lines, 0((\eB\/m 2 ) n ), thereby 
one has to sum up all the diagrams when B is larger than the critical field B ~?t B c = 
m 2 /\e\ to obtain a "dressed" fermion propagator. This effect is shown in Fig. [TJ Since this 
dressed fermion propagator includes all-order contributions with respect to eB, any process 
involving the dressed fermions becomes nonlinear with respect to the external field B, and 
nonperturbative. Similarly, one can define the critical electric field E c = m 2 /\e\ which also 
indicates breakdown of the ordinary perturbative expansion. However, electric fields beyond 
E c induce instability of the vacuum j^, [2^, called the Schwinger mechanism, and will be 
screened by creation of fermion-antifermion pairs. On the other hand, the critical magnetic 
field B c just indicates onset of strong nonlinear effects, and it makes sense to treat magnetic 
fields stronger than the critical one. Thus, as far as one considers static magnetic fields, 
one can discuss very strong nonlinear regime in QED. This is called the "nonlinear QED" 
regime. 



If there are several species of fermions iftw with different masses mm and charges e^, 
one can define critical magnetic fields = m f^/\^(i)\ corresponding to each fermion: 
BP < BP <■■■. Since magnetic fields can become stronger than the (minimum) critical 
field (unlike the critical electric field beyond which the Schwinger mechanism occurs), one 
has to include all the relevant fermion degrees of freedom when one treats very strong 
magnetic field; typically if Be < B < B^ 1 ^, then at least we need to include fermions up 
to ip( n \ In the present paper, however, we provide formula only with the lightest fermion. 
Contributions of heavier fermions are easily obtained by appropriate replacements of the 
electric charge and the mass by those of the heavier fermions. Note also that contributions 
from different fermions are all additive in the vacuum polarization tensor. 

A fundamental key ingredient is the fermion propagator in an external field G(x, y\A c \) = 
(0\Tijj(x)ijj(y)\0) which includes all the interactions with the external field. Recall that the 
propagator is a Green's function of the Dirac operator. If one includes only the external 
field into the Dirac operator, the propagator satisfies (in the momentum space), 

(f - e4d ~ m) G(p\A c} ) = i . (2) 

Therefore, one finds 

G{p\A d ) = ' (3) 

p — e4 c i ~~ m 



oo 

% 



n=0 



rv 



(4) 



The last expression allows us to diagrammatically depict the nonlinear interaction with the 
external field as shown in Fig. [TJ 

When one considers the propagation of a photon in a strong background field A^, one uses 
the fermion propagator ([3]). For example, a one-loop contribution of a photon self-energy 
is computed as in Fig. [2j Thus, modification of the fermion propagator due to the external 
field causes a significant influence on the photon propagation. Representing the photon 
self-energy, or equivalently the vacuum polarization tensor, in a constant external field as 
i x ~ u)i one finds that the kinetic term of the propagating mode a M acquires a quantum 
effect as, =Sf\ oop = | fdy a M (x) n^(x — y) a v {y). Then, the Euler-Lagrange equation for 
a M leads to a modified Maxwell equation, [ c^rf v — q fi q u — Il^'(g) ] a u (q) = 0. Although the 
vacuum polarization in the ordinary vacuum does not modify photon propagation owing 



FIG. 1. Dressed fermion propagator (a double FIG. 2. One-loop diagram of the 
line) includes all the tree-level interactions with a vacuum polarization tensor in a 
strong external field (wavy lines). strong magnetic field. 

to the gauge and Lorentz symmetries, we show that quantum excitations in an externally 
applied electromagnetic field behave like electron-hole excitations in dielectric substances 
3] . Indeed, the vacuum polarization tensor corresponds to a response function of the Dirac 
sea to an electromagnetic field induced by an incident photon. 

To compute the vacuum polarization tensor in the external field, we use the "proper-time 
method" which was developed by J. Schwinger 24J j . One can equivalently rewrite the dressed 
propagator ([3]) in a different way as 

G(p\A cl ) = i{f- e4d + m) x - / dr j*{&-*4*) a -{m*-*)} ? ^ 

1 Jo 

where the integral with respect to the "proper time" f is convergent owing to a prescription 
by —%£. Note that f has dimension of inverse mass squared. Infinite sum with respect to the 
external field is now encoded into the exponential factor in the integrand. As summarized 
in Appendix El one can explicitly compute G(p\A c {) when the external field is constant. 

By using the dressed propagator shown in Eq. (1A7|) . we can now calculate the 1-loop 
vacuum polarization tensor in the external magnetic field (see Fig. [2]): 



i IC(g) = He) 2 (-1) 



d 4 p 



Tr ^G(p\A cl )YG(p + q\A c 



(6) 



where an overall minus sign arises from the fermion loop. This is a gauge invariant quantity 
and should be independent of the gauge we adopt Note that momentum integration gives 



rise to an ultraviolet divergence, because the fermion propagator behaves as G(p\A, 



V 



for large p. While we find a quadratic superficial degree of divergence from a naive power 
2 In the actual calculation, however, we adopt specific gauges. For computation of the fermion propagator 
G(p\A c \), we worked in the Fock-Schwinger gauge for the background field A c \ (see Appendix |A")) . As we 
will see below, we then compute the polarization tensor in the covariant gauge which fixes the residual 
gauge symmetry of the dynamical gauge field a M . 
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counting, it actually diverges logarithmically owing to a specific tensor structure demanded 
by the Ward identity, q^H^^q) = 0. We will explicitly show a gauge invariant form of the 
vacuum polarization tensor in the following, and comment on renormalization in the end of 
this section. 

It has been known that the vacuum polarization tensor is endowed with a gauge-invariant 
tensor structure in terms of three transverse-projection operators satisfying q^P^ = 0. Their 
explicit forms are given by 

iC(?) = - (xo*r + XiPT + X*PT) , (7) 
PZ" = <?rr-M, iT = - W ' = ( 8 ) 

where we suppressed arguments of Lorentz-scalar functions Xi (* = 0, 1, 2), of which explicit 
expressions are immediately shown below. We now apply the external magnetic field along 
the third axis of spatial coordinates in the negative direction so that eB 3 > for an electron 
(B 3 is the third component of the magnetic field vector B l ). Since the presence of the 
magnetic field specifies a preferred direction, one distinguishes longitudinal and transverse 
directions with respect to the magnetic field. The metric tensor if = diag(l, —1, — 1, —1) 
is decomposed into longitudinal and transverse subspaces rj^ = diag(l,0,0, —1) and 77^ = 
diag(0, —1,-1,0), respectively. Similarly, longitudinal and transverse momenta are defined 
as = (g°,0,0,g 3 ) and q± = (0, 9 1 , g 2 , 0), respectively. 

In the proper-time method (see Appendix|X]for more details), the polarization tensor (JH]) 
has integrals with respect to the momentum p, and two proper times which we denote t\ and 
f 2 corresponding to each propagator (see Eq. (JSj)). The integration over p is just a Gaussian 
integral, and can be straightforwardly performed. Then, we rewrite the remaining integrals 
by using two dimensionless variables r = eB[j\ + Tz)/2 and (3 = eBiji — t 2 )/t. Notice that 
the scalar functions Xi (i = 0? 1)2) are dimensionless, and thus we further introduce three 
dimensionless variables, B v = B/B c , rjj = and r\ = = — so that the scalar 
functions Xi are expressed as 

X t (rl rl; B t ) = ^ ! d/3 / °dr e - i?tcos(/3r) e^ cotT e~^ T , (9) 

47r J-l Jo SUIT 

where we have introduced two shorthand notations: rj = —2 r\ IB T and u = 77 /sin r, and 



lastly 0|i and Tj are known functions given by the following 



14H17J, 



19| (see also Refs. 



20 



for details): 



rf,5 r ) = -Ul-(l-/? 2 )rf}, (10) 



and 

r (r, f3) = cos(/3r) — (3 sin(/3r) cot r , 

r 1 (r,/3) = (l-/3 2 )cosr-r (r,/3), (11) 

cos(/3t) cost 

J-2(t,/}) = 2 — i-o{r,p) ■ 

sin r 

In Eq. (Q, the coupling constant in the overall factor a = e 2 /4ir comes from propagating 
photons attached to the fermion one-loop, while the others are from higher-order effects 
associated with the external magnetic field as is evident from the observation that they 
always appear with the magnetic field as eB. As will be shown in Sec. IIVI the scalar function 
Xi is a real- valued function when r 2 < 1, reflecting convergence of the double integral in 
Eq. (Q, and the same holds for the others, Xo an d X2> when r 2 < + + 2B r ) 2 J A. Whereas 



reliable numerical computation has been performed in this kinematical region |18|, analytic 
calculation is necessary for understanding behavior of the polarization tensor in the whole 
kinematical region. In Sec. IIVI we will explicitly perform the remaining two integrals to 
obtain the analytic expression of the polarization tensor. 

We emphasize here that we have neither specified any dispersion relation for the external 
photon momentum, nor even assumed that the photon is on-shell. We will discuss later 
that the photon dispersion will be determined as a result of interaction with the external 
magnetic fields. In general, there will be a deviation from the massless-type dispersion 
relation indicating that rjj+r^ ^ and u ^ \q\. Accordingly, the longitudinal and transverse 
momenta appearing in Eqs. (J9J) - ( ITTj) should be treated independently at this moment. 

Recall that only xo survives in the ordinary vacuum and is divergent. Divergence of Xo 
is also seen in the presence of the magnetic field (while the others are finite). This can be 
explicitly verified as follows. First of all, let us see the lower limit of r integration in Eq. (j9]). 
If one takes the limit r — > in Eq. (TTTT) . one finds two of them vanish as T\, I^ — > and thus 
Xi and X2 are finite, while xo is logarithmically divergent because ro — > 1 — /3 2 as r — > 0. In 
other words, subtraction of To(t, (5) in the definition of r\ and I^ ensures the finiteness of xi 
and X2- Next, one can easily identify the origin of the logarithmic divergence by referring the 
diagrams in Figs. [1] and [21 If one resolves the dressed fermion propagator (a double line) in 

8 



Fig. [2] into the sum of external magnetic field insertions as shown in Fig. (TJ the first diagram 
is given by the same vacuum polarization tensor as in the ordinary vacuum without any 
external field attached. This diagram has a logarithmic divergence as conventionally known. 
We also notice that this is the only divergent diagram in the series because any additional 
external leg reduces the degree of divergence. Therefore, the vacuum polarization tensor 
(J9]) contains the same divergence as in the ordinary vacuum, without any additional one. 
As shown in Appendix [B] in detail, we can remove this logarithmic divergence by adopting 
an on-shell (with zero transverse momentum) renormalization conditional such that a finite 
vacuum polarization tensor tl^(q\\, q±; B r ) vanishes at the simultaneous limits of on-shell 
and vanishing field, r 2 ,r\ = and B T = 0. On the basis of this prescription, we define a 
finite vacuum polarization tensor, 

U^(q h q±; B r ) = IL£(q hq±] B r ) - 11^(0) , (12) 

where the vacuum polarization tensor in the ordinary vacuum n.Q U (q 2 ) is obtained in the 
vanishing field limit as liniB^o IT^(<?||, <?±; -Br) — > n(T(?)- We will more specifically discuss 
these things related to renormalization in Appendix [B] 



III. OPTICAL PROPERTY OF THE VACUUM WITH MAGNETIC FIELDS 

Before we go into the detailed calculation of Xii let us discuss observable effects of nonzero 
Xii which does not require explicit expressions of Xi- I n particular, we introduce the notions 
of dielectric constants and refractive indices that are frequently used in context of optics of 
dielectric substances and semiconductors. 

Propagating photons in the ordinary vacuum have two transverse oscillating modes with 
the dispersion relation of the massless type, u 2 = \q\ 2 , which does not vary even if quantum 
corrections are included, as long as the gauge and Lorentz symmetries are preserved. The 
presence of an externally applied magnetic field however breaks the Lorentz symmetry, and 
the dispersion relation of a photon is subject to modification which leads to two intriguing 
phenomena: Two propagating (physical) modes can have refractive indices different from 
unity, and a photon is able to decay into a fermion-antifermion pair if its energy is large 

3 Note that the divergence has not been regularized in Eq. (j9]) - (THJ. We will introduce appropriate 
regularizations depending on analytic and numerical calculations, which will be shown in Appendix [Bl 
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enough. In this section, we show kinematical aspects of these effects on the basis of the 
Maxwell equation modified by a nonlinear interaction with the external magnetic field. 

In the presence of an external magnetic field, the Maxwell equation for the dynamical 
photon field a M is modified as 

1 



q 2 rf v 



gV-nSTfai.ffjL;^) 



a v (q) = 



where £ 9 is a parameter in the gauge fixing term, JzfcF = — ^-(d^a^ 
vacuum polarization tensor ([7]), one finds 



a v (q) = 



(13) 

Substituting the 
(14) 



Below we explain in a qualitative way how the vacuum birefringence and on-shell photon 
decay follow from this modified Maxwell equation. To this end, explicit expressions of the 
scalar functions Xi are n °t important. 

In order to identify two physical modes in the modified Maxwell equation ( fl4|) . we need 
to 'diagonalize' the matrix equation by selecting appropriate vectors associated with the 
direction of the magnetic field. First of all, it should be noticed that one can define three 
independent vectors that are orthogonal to the momentum of the propagating photon g M . In 
addition to the invariant tensors r/^ and ^ vpa , one can utilize the field-strength tensor F^ v 
and the photon momentum g M to construct independent vectors. Namely, one finds that the 
following four vectors (A = 0, ■ • • ,3) are orthogonal among each other [251 ] : 



(o) 
i-i 

(i) 



q 2 F^F up qP - {q v F vp F<>° q°)q» 



F^q 



V {2) 



F» v q v 



J (3) 



(15) 



where the dual field-strength tensor is defined as F pu = \e lll ' p<T F pa with the completely 
antisymmetric tensor, e 0123 = 1. Indeed, one can show the orthogonal relations 3: «a) v {a)n = 
for any A and cr(^ A). 

Assuming that there is only a magnetic field and it is oriented to the third direction of 
the spatial coordinates, the set of vectors f|T5]) simplifies to 

4 In general, when electric (E) and magntic (B) fields are present, one finds vf^ v (ij)Lt = except for 



(A) »^ 

= —q 2 (E ■ B). Therefore, the orthogonality among all the vectors holds when E ■ B = and 



thus E = (our case). 
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vfc = Bft, < 2) = 5gj\ vfc = <f, (16) 

where B is a magnitude of the external magnetic field. For convenience, we have introduced 
two momentum vectors qi^ = (q 3 , 0, 0, q°) and q^_ = (0, q 2 , —q 1 , 0). They are orthogonal to all 
the previously introduced momenta, g M , q^ and q±. Therefore, the vectors and q± span a 
space complementary to the space spanned by qff and q^_. 

Notice that the projection operators in the modified Maxwell equation (j!4p are orthogonal 
to the momentum vectors, q^, q^ and q±, namely, q^P^ = q\\^P i fJ,u = q^P^ = 0. This 
fact suggests that one can represent the projection operators in a different way by using the 
vectors q^ and q^_. Indeed, two of the projection operators are expressed as P-^ v = —q^q^ 
and P 2 ^ v = q±q±- Therefore, one finally finds that the projection operators and thus the 
modified maxwell equation are 'diagonalized' as follows: 

?r = i 2 Ko)<o) + <d^(i) + 4)^(2)) , pr = 4 > p ^ u = £ > 

3 

M^{q)a v = ^M(A)(g)7rf A) 7r^ a v = , (18) 

A=0 

where we have rescaled Vq\ so that they have unit norms ajf| 



(A) 



(A = 0,1, 2, 3). (19) 



and the coefficients M.(\) are expressed in terms of Xi- One can easily find the eigen- modes 
of the modified Maxwell equation by expanding the dynamical photon field in terms of ir^ 
and assuming a plane wave solution: 

3 

a»(q)=N < A) e-^-^ , (20) 

A=0 

where q^ = (u, q) and N is a dimensioned normalization constant. Notice again that we 
have not specified the dispersion relation which provides a relation between u and q, and 
this is what we suppose to obtain by solving the modified Maxwell equation. Indeed, one 



Our definition of ir^ contains somewhat confusing notation wherein the subscripts of the last two projec- 
tion operators in Eq. (|17l) oppositely correspond to those of the vectors on the right-hand sides. However, 
we keep this notation because the subscripts originally defined in Eq. (|15p give natural ordering (i.e., Pi 
is given by ftm, etc) in case of electric fields. 
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finds that linear independence of ir£<. leads to a set of eigenvalue-equations 



{(1 + Xo) q 2 }^ 0) =0, 

{(i + Xo)g 2 + x 2 gi }^ = o, 

{ (l + Xo) ? 2 + Xi gj } 7rf 2) = 0, 



which immediately yields the dispersion relation of a photon: 

\q\ 2 



CO 2 



(22) 



Here we have defined a dielectric constant e which persists unity in the ordinary vacuum. 

It is easily seen that the first and fourth equations lead to dispersion relations of the 
massless type, u 2 = |q| 2 , with the speed of light in the ordinary vacuum (e = 1), as long as 
neither (1 + xo) = nor ^g 1 = 0. However, as explained in Appendix [Dj these two modes 
turn out to be unphysical (see Eqs. (ID3P and (lD4j) ) and thus we do not discuss them below. 

In contrast, the second and third equations yield nontrivial dispersion relations of physical 
modes. Without loss of generality, we suppose that a photon is propagating in a plane 
(y = 0) spanned by the first and third directions of spatial coordinates, and that the external 
magnetic field is oriented to the (negative) third direction. The momentum vector is now 
represented as = (u, \q\ sin(7r — 6), 0, \q\ cos(7r — 6)), where 9 denotes an angle between the 
direction of the external magnetic field and the momentum of a propagating photon. Then, 
the second and third equations yield the vacuum dielectric constants e± and en, respectively: 

e _L = ~ — ^° fc> r TTm j (23) 

1 + X0 + Xl for < 2V (24) 



1 + Xo + Xi cos 



4 W 



19] in the radiation gauge, while we have 



The same expressions were obtained in Ref. 
obtained them in the covariant gauge. The coincidence of two calculations in two different 
gauges indicates gauge invariance of the dielectric constant, as expected. Whereas we have 
not specified the explicit form of Xi, two constants e_|_ and e» are in general not equal to unity 
and different from each other[] Namely, we have found that the externally applied magnetic 



6 We note that the same results can be obtained from the pole position of a photon propagator in the 

magnetic field (sec Appendix [C]) . 
' In the ordinary vacuum, only \o is nonzero. Even with nonzero xo, the above expression reduces to 

e±_ = en = 1 when \i = X2 = 0, as expected. 
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field gives rise to two distinct dielectric constants depending on the polarization directions. 
This immediately implies that we can define two distinct refractive indices. Therefore, we call 
this phenomenon vacuum birefringence after a similar phenomenon in dielectric substances. 
Notice also that, due to the violation of the Lorentz symmetry by an external magnetic field, 
the dielectric constants explicitly depend on the zenith angle with respect to the magnetic 
field direction. Nevertheless, the system maintains a boost invariance in the direction of the 
constant external magnetic field, and thus photon propagations in the directions at 9 = 0, ir 
are special. Substituting these angles into Eqs. ( )23|) and ( 124)) . we find that both of the 
dielectric constants become unity e±(9 — 0,7r) = en(6 = 0,7r) = 1, as a consequence of the 
boost invariance. 

Let us further assume that the scalar functions x% have imaginary parts. Then, the 
dielectric constants also have imaginary parts: e = e roa i + i eimag- We will see later in Sec. IIVI 
that imaginary parts indeed appear under some kinematical conditions. Since the dielectric 
constants and refractive indices are related with each other via n 2 = e, we can similarly 
define real and imaginary parts of the refractive indicesjf 



n 



n vc&\ 4" i "imag j (25) 



where one can easily find (|e| = Je 2 cal + ef ) 



™real = ~^V\ 6 \ + e real, (26) 
1 



^irnag = l e l ~~ 6 rcal • (27) 

While a purely real dielectric constant yields a real refractive indeXG ^reai = y/^reai, a purely 
imaginary dielectric constant leads to an equal magnitude of real and imaginary parts, 
n rC ai = '"■imag = \/l e l/2 with |e| = |e imag |. In Appendix [Dj we have explained that real and 
imaginary parts of the two refractive indices nu and n± in general provide distinct phase 
velocities and extinction coefficients of physical modes A = 1,2 in Eq. (120]) . respectively. 

It should be noticed and taken care of carefully that the right-hand sides of Eqs. ( )23l) 
and (1241) depend on the dielectric constants e± and e\\, respectively, through the photon 

8 Since there are two polarization modes, one can define two different refractive indices n\\ and n±, both of 

which can be complex quantities. 

9 Inserting Eqs. ([23]) and (|24[) into this expression and then expanding them with respect to xo ~ 0, 

Xx "C 1, and X2 < 1, we obtain approximate forms shown in literatures, n± = 1 — (x2/2)sin 2 and 



1 = 1 + (xi/2) sin 9. However, one should bear in mind that xo, Xi an d Xi can t> e even divergent and 



complex depending on kinematical variables, as w iU see m Sec. IIVI 



momenta rjj and r\. According to the definition of the dielectric constant (|22|) . those 
momenta are written as r 2 = u) 2 (l — e cos 2 9) and r\ = — e a) 2 sin 2 6, where we introduced 
a normalized photon energy Co 2 = u 2 /{Am 2 ). Thus, if one wants to obtain the dielectric 
constant e as a function of the photon energy u), Eqs. f l2Hj) and f[2"4"j) have to be solved 
with respect to e. Physically, this structure reflects effects of a back-reaction appearing as 
screening of an incident photon field by an induced vacuum polarization, which should be 
taken into account self-consistently when the magnitude of the polarization becomes large. 
As mentioned above, the dielectric constants are in general complex, and thus damping of 
the incident photon field, due to the decay into a fermion-antifermion pair, should also be 
treated self-consistently. Therefore, we need to simultaneously solve two sets of coupled 
equations (obtained from Eqs. ( I2"5j) and fl2~4"j) ) for the real and imaginary parts of e± and ey 
in a self-consistent way. There are indeed kinematical regions where these procedures play 
an important role to correctly obtain the dielectric constant. We will explicitly demonstrate 
it in the next paper 22]. 



IV. ANALYTIC EVALUATION OF THE VACUUM POLARIZATION TENSOR 

In the previous sections, we have discussed in a formal way that a nonzero vacuum 
polarization tensor in a magnetic field necessarily gives rise to the vacuum birefringence. 
In order to know how large the effect is and how it depends on kinematical conditions, we 
have to explicitly evaluate the scalar functions Xi (i = 0, 1,2) of which formal expressions 
were already given in Eqs. fl9]) - (II ip . Notice that this representation contains somewhat 
complicated integration with respect to (3 and r, corresponding to the difference and the 
average of two proper-time variables. In fact, this complexity has prevented previous studies 
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191 ] . and even 
We 



from complete analytical understanding of the vacuum birefringence 
from performing numerical computation except in a limited kinematical region 
will, in this section, perform the integration analytically by rewriting the integrand into 
a double infinite series of familiar functions, and show that the infinite series indeed has a 
physical meaning. We will see that the scalar functions x% an d thus the refractive indices are 
sensitively affected by microscopic structure of the fermion spectrum in a magnetic field, i.e., 
the Landau levels formed by the fermions. Therefore, as in optics for materials, macroscopic 
electromagnetism of the vacuum should be elucidated by the microscopic dynamics of the 
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Dirac sea which is excited by an incident photon field. 



A. Computing the scalar functions \i 

Before we present the procedures how to analytically perform the integrals, let us briefly 
point out technical difficulties in the integral representation, Eqs. fl9]) - fTTTj) . Note that each 
integrand in \% contains an exponential function whose argument is given by trigonometric 
functions. This factor causes strong oscillation of the integrand composed of arbitrarily 
higher harmonics, because k-th power of the trigonometric function that appears in a Taylor 
expansion of the exponential factor is composed of higher harmonics up to k times the 
fundamental frequency. For example, the factor e J??cotT contains a term proportional to cos 3 r 
at k = 3 in the Taylor expansion, which is composed of up to the 3rd higher harmonics as 
cos 3 r = {3 cost + cos(3t)}/4. In general, this strong oscillation even cannot be periodic 
due to mixing of two fundamental periods, 2ir and 2-k//3. These complicacies, directly or 
indirectly, cause difficulties in analytic and numerical evaluations of the integral. In fact, 
as we will see later, results of the integration have singular behaviors at some kinematical 
conditions, which may invalidate numerical approaches. 

The first thing to overcome this situation is to rewrite the integrand so that the exponen- 
tial function does not contain trigonometric functions of the fundamental period, 2tt/(3 (Step 
I). We will still have an exponential function with a trigonometric function in the argument 
which is however composed of a unique fundamental period, 2ir. We can rewrite it by a se- 
ries expansion so that the exponential function does not contain any trigonometric function 
(Step II). The resultant integrand will allow us to easily perform analytic integration (Step 
III). Below, we explain these procedures step by step. 



1. Step I: Partial wave decomposition 

Consider the functions in the integrand in Eq. ([9]) that depend on trigonometric functions 
with a fundamental period 2n/f3. Namely, e - iMCOS ^ T ), cos(/3r)e- incos(/?T) , and cos(/3r)e- mcos(/?T) 



We can equivalently rewrite these functions by using the partial wave decomposition 



26). 
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By adopting the formulae given in Appendix [Ej we obtain the following^: 

pi poo 00 o\ 

Xl = il / d/3 / dT Y (2 - £ n0 ) * 1 e ^cotr e -i(0 |r n^)r ^ 

4 W-i io smr 
with 7,- n ' ) (r, /3) defined by 

7^ n) (r, /3) = -{J n+ i(-m) + I n _ x (-iu) } - n/3 r]~ l I n {-iu) cost, 

7 i" } (r, /3) = (1 - /3 2 )4(-^) cos r - ^ n \r, (3) , (29) 

72™ } (r, /3) = sin~ 2 r { I n+l (-iu) + J„_i(-m) - 2I n (-iu) cost } - ^(t, (3) , 

where I n (—iu) is the modified Bessel function of the first kind. Except for the last trivial 
exponential factor e~ l ('l > \\~ n ^ T in Eq. fl28|) . the integrand now oscillates with the fundamental 
period 2ir. Note that I n (—iu) also oscillates with the fundamental period 2ir because its 
argument is proportional to cosecant as, — iu = —ir/cscr. 



2. Step II: Infinite series expansion 

We are able to further simplify the integrand by using the following technique. We first 
introduce a variable z defined by 

z = exp(— 2ir) . 

Then, a product of the exponentiated cotangent and the modified Bessel function in Eqs. ( 1281) 
and ([29]) are rewritten as e ivcotT I n {-iu) = e~ v e~ z ^ J n ( 2( f_y ). We use the fol 



formula which is deduced from the relation between confluent geometrical functions 



owing 



exp 



3: 

4 i^f) =&-.) (*»*)* t * , (30, 

where L™(x) and T(x) are the associated Laguerre polynomial and the Gamma function, re- 
spectively. Notice that this relation disentangles the function with a nontrivial z-dependence 
(left-hand side) into a sum of polynomials of z (right-hand side), which is quite useful for 
the present purpose. Namely, the resulting integrand can be now written only in terms 
of simple exponentials of r. Taking parameters in the above formula as x = y = t] and 



10 We just use the partial wave decomposition as a mathematical tool. Later, the index n is interpreted as 
a physical quantity, but not an angular momentum. 



1(3 



using (1 — z) 2 = YliLo ^>T=o = ^fc=o(^ + 1) z * C ' we ^ nc ^ that ^he double integrals are 
simplified to either of the following representations: 

F?{rf u B T ) = 4- /°dr e - < (*n+ 2< -^+ B ) r , (31) 

B r J -I Jo 

G?(rjf,£ r ) = -^- / dp [°° dr Pe- i fa +2t - nf>+n ) T , (32) 

B r J -I Jo 

/l /*oc 
d/3 / rfT / 52 e - l (^ || +2£-n/3+n)r > (33) 
-1 Jo 

Appart from these integrals, we have nontrivial coefficients which contain the associated 
Laguerre polynomials as we will explicitly show in the final expression below. Here we just 
emphasize that, thanks to the formula (130]) . the longitudinal (rj?) dependence of Xi( r y; r i) 
appears only in the functions F, G and H, while the transverse (rj_) dependence are factorized 
and appears in the other parts. 



3. Step III: Performing double integrals 



The integrals in Eqs. ( 13111 - ( 1331) are simple enough and we can carry out integration 
analytically. However, it should be noticed that these integrals contain singular behaviors, 
which appear when the exponential factor becomes unity due to a vanishing argument. 
Depending on the photon momentum r 2 , there arises such a divergent contribution within 
the integral region, /3 G [—1,1]. A simple way to obtain the integrals with respect to any r| is 
provided by the technique of analytic continuation. First, assuming that the rescaled photon 
momentum r| satisfies <j)\\ (r 2 , B T ) + (2£ — n(3 + n) > 0, we perform convergent integration by 
rotating the contour in the complex r-plane downward. Then, analytic continuation allows 
us to inspect the analytic property afterward. At this moment, the r- integral is carried out 
as follows: 



i 



dp 



m^Bj-J^ r p2_ nB ^ +{1 _ r 2 ) + {2e+n)Bi = J M(i) • (34) 

The integral with respect to P simply provides a difference I^{rT) = I™ + (r 2 ) — /™_(r|) 
between the following functions: 

/; ( r 2) = ^ =arcta n ( _A ±J^_ ) . (35) 
t±K V y/Aac - b 2 W4ac-by V ; 

We have introduced parameters a = r 2 , b = —nB r and c = (1 — rj?) + (2£ + n)B T , and the 
integral with respect to P has been carried out in a regime, 4ac — b 2 > 0. Note that the 
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condition for r 2 assumed above involves this regime. Similarly, the remaining other two 
integrals are performed as IH. 



G?(rjj,5 r ) 



(36) 



nB T 



^(Br) + — 2 { (b 2 - Aac) + {nB v f } I t 
4r ll 



£AK r \\ 



(37) 



where we have introduced a real constant 5 (i.e., independent of kinematical variables): 



3?(B r ) = In 



1 + 2£B r 



1 + 2(£ + n)B r 



In 



m 2 + 



2(£ + n)eS 



We have performed the integrals under the condition 4ac — b 2 > 0. However, the obtained 
results can be analytically continued to the regime 4ac — b 2 < 0, where the function -/^(r 2 ) 



will become complex as we will see later. 



4- Final results 



Now that we have performed all the integrals, we find the analytic expression of x% 
(i = 0, 1, 2). In order to present the results in a compact form, let us define the coefficients 
C™{rj) expressed by the associated Lagguerre polynomial^ : 



(38) 



(£ + n)\ 

Then, the coefficients Xi of the vacuum polarization tensor (JTj) are finally represented as the 
infinite sum of known functions as 

aB T 



Xi 



Air 



E ( 2 - 6 ^ 



n=0 



n(2) 



1=0 



1=2 



(39) 



where we have introduced three functions and Q™f\ as we define below. 



For % = 0, they are defined as 



n(0) 
£0 

n(l) 
^0 



;i - s n0 ) cr\v) FP(^ B T ) - nrj' 1 C^v) B t ) , 



11 Although the expressions of Gr"(r 2 , B r ) and U"(r 2 , B r ) apparently look singular at r?, = 0, one can check 



I 

that both of them arc regular there. 

12 We employ a convention for the associated Laguerre polynomial defined by series: L"(rj) 

i 

r=0 
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= o . (40) 

Note that the functions F and G (and also H) depend on r| only through the combination, 
£ = 2r?JB T . Thus, we have written them as F(£,B T ) and G(£,B T ). Similarly, for % = 1, we 
have 

^ 0) = cnv) wee, s r ) - ^)) - < 0) > 

fi;i 2) = o, (4i) 

and finally for z = 2, 

Q n(0) _ _ o n(0) 
J£ £2 — , 

^ 2) = D^ 2 \rf) Fj?(£, £? r ) . (42) 

In the last expression, we introduced coefficient functions .D 1 - 1 ) and defined by 

e-i 

= -8 E (^-^{(l-^GT 1 ^)"^)} , (43) 

A=0 
1-2 

= -»£ (^-A-ijia+ycf 1 ^)-^)} . (44) 

A=0 

These are the final analytic results represented as an infinite series of known functions. For 
fixed n and £, each term in the expansion is finite (except for some special kinematical points 
as we will discuss soon). The ultraviolet divergence of the function xo appears only after 
one takes the infinite summation (see Appendix iBl) . 

Figure [3] shows ^-dependence of the coefficient C™(rj) in Eq. (155]) at fixed rj and n. The 
associated Laguerre polynomial provides an oscillation, and, in case of finite n, a ratio of the 
factorials strongly suppresses C™{rj) for small £. In Appendix [F] we discuss the behavior of 
C™{rf) in the limit of large £ ^> 1. By using the limiting form (1F5I) of Cf(rj) at large values 
of £ with fixed r\ and n, one can put the upper bound for C™{rf): 



< Cyfa) < — ^= e" 1 ^ (I > 1) , (45) 

Black dashed line shows the upper bound of C"(r/) for n — 0, and the other two cases are 
also bounded at large £ outside the plot range. While the upper bound in Eq. (J4"5j) indicates 

19 




FIG. 3. Coefficient C™(r/) at fixed n and r\: C™(rj) at discrete £ takes a value on curves plotted 
against continuous I. A dashed line shows the upper bound for the coefficient at n = that is 
valid at large I 3> 1. Owing to Eq. (|F5p . other two curves are also bounded at large I outside the 
plot range. Note that the coefficient has a finite value at £ = given by Cq (rj) = e -T? 77 n /n!. 

that C™(rj) for large I converges to zero by an inverse-square- root, this suppression is not 
sufficiently strong to guarantee a convergence of the infinite sum Indeed, we find a 

logarithmic divergence of the function xo by carrying out the infinite sum. However, the 
logarithmic divergence has been foreseen a priori to appear in terms of the ultraviolet diver- 
gence mentioned below Eq. ( fPT|) . which should be removed by an appropriate prescription 
(see Appendix iBl). 

B. Physical meaning of the results: Landau levels 

The use of infinite series expansion was motivated rather by technical reasons for making 
the double integrals simpler. However, as we will see below, the indices i and n in the 
expansions have well-defined physical meaning: they are related to the Landau levels of the 
fermion-antifermion pair which appears in the one-loop diagram shown in Fig. [2j 

First of all, it should be noticed that an important kinematical property of the scalar 
function x%> with respect to rj?, is essentially determined by the function /^(rj?) which is 
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the only one possible source of an imaginary part. Therefore, we investigate the structure 
of 

-^A( r ii) m or der to understand kinematics in the presence of an external magnetic field, 
and to identify the physical meaning of the indices i and n. 

A kinematical property of is specified by a discriminant appearing in Eq (I35j) : 



V = b 2 - Aac 

= (-nBr) 2 - Ar\ [(1 - rf) + (2£ + n)B t ] 



(46) 



Depending on the sign of T>, arguments of the arctangent in If^{rT) are either real or pure 
imaginary. Solving an equality, T> = 0, we find boundaries of the momentum region specified 
by a pairwise solution, 



- { v 7 ! + 2£B r ± y/1 + 2{t + n)B r y = s% 



(47) 



If the photon momentum resides in the regime sfT < r 2 < s+ where T> > 0, the function 
7^(r|) is obviously a real- valued function. On the other hand, in the other regimes, rjj < s e ™ 
and sfp < r 2 where T> < 0, it possibly becomes a complex function. However, we show in 
Appendix [G] that If^{rT) has an imaginary part only if the photon momentum resides in 
the higher regime, s+ < r 2 . In other words, it does not have an imaginary part in the 
complementary regime, r| < s!p which includes rjj < s e ". After careful investigation shown 
in Appendix [Gj we obtain a piecewise representation as, 



i 



/( r 2_ s ^)( r 2_ s fn) 
1 



| In 



a— c— Vb 2 — iac 



a—c+\/b 2 —4ac 



in 



^)( r 2_ s <jn)| 
1 



arctan 



b+2a 



'W-si n )M-sV 1 ) 



111 



\/Aac—b 2 
a~c—\/b 2 —4ac 



a—c+\/b 2 —Aac 



arctan 



2m 



b-2a 
\jAac-b 2 



[rf, < s!T 



(si" < r 2 < s<£] 



(4? 



(sfp < rj?) . 



An imaginary part in I^{rT) results in the imaginary part of the scalar functions \ii 
and of the vacuum polarization tensor U^(q). Consequently, the dielectric constants and 
thus refractive indices have imaginary parts as alluded in Eqs. (1231) - (127]) . They represent 
damping of photon propagation in the presence of an external magnetic field. As foreseen 
owing to the optical theorem, this damping is caused by decay into a fermion-antifermion 
pair. Thus, we interpret r 2 = s+ as the threshold momentum of photon decay, and notice 
that there are infinite number of threshold momenta as the indices I and n run from zero to 
infinity. 
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In terms of dimensionful quantities, the threshold condition rjj = s+ for the photon decay 
is expressed as 

g j = | Vm 2 + 2£eB + ^m 2 + 2{£ + n)eB | 2 . (49) 

Recall that the energy of a charged particle in a magnetic field is £ n (p 2 ) = \/m? + p z + 2neB. 
This immediately implies that the right-hand side of Eq. (1491) exactly agrees with the in- 
variant mass of a fermion-antifermion pair carrying quantized transverse momenta in the 
magnetic field and vanishing longitudinal momenta along the external field. The integers £ 
and £ + n specify the Landau levels. 

Reflecting the boost invariance of the system along the external constant magnetic field, 
the left-hand side of Eq. (1491) has a boost invariant form g 2 . It indicates that the decay 
condition should be invariant in a class of Lorentz frames which are connected by the boost 
along the external field, because the propagating photon receives an influence from the 
magnetic field of the same configuration. Owing to the boost invariance, it is possible, 
without varying the configuration of the external field, to take the Lorentz frame in which 
the longitudinal momentum of the photon vanishes! 13 ! q^ z = 0. In this Lorentz frame, we 
have g 2 = — q^ z = w 2 ^, and the condition (T49J) represents nothing but the smallest 
photon energy to produce a fermion-antifermion pair in the Landau levels £ and £ + n with 
vanishing longitudinal momenta. 

For each £ and n, the function /™(?" 2 ) has a single threshold momentum beyond which 
a photon can decay into a fermion-antifermion pair with the Landau levels specified by £ 
and £ + n. Since the coefficients Xi an d thus the vacuum polarization tensor are given as 
the infinite sum of /"(r 2 ) over the indices £ and n, there are infinite number of threshold 
momenta in the photon momentum r 2 . Beyond the each threshold momentum, the vacuum 
polarization tensor has an imaginary part, indicating a branch cut continuously running on 
the real axis in the complex gjj-plane. The half line of the cut and its starting point at the 
threshold correspond to the continuous and vanishing longitudinal momenta of a real pair 
excitation from the vacuum, respectively. 

In order to confirm our interpretation of the Landau levels, let us look at the contribution 
from the lowest Landau levels, £ = n = 0. Taking £ = n = in Eq. ( 1391) . we find that the 



13 If one wants to perform an equivalent analysis in a Lorentz frame disconnected to this class, one has to 
take into account an electric field orthogonal to the magnetic field. 
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coefficients Xi °f the vacuum polarization tensor survive only for 2 = 1. Namely, they read 
Xo = X2 = 0, and 



2 ri;S r ) 



r e-xi,{/0 A (rJ)-2} 



47T 



(50) 



where we have used L^rj) = 1. Depending on the values of r 2 (the boundaries are s_ = 
and = 1 (or, g 2 = 4m 2 ), see Eq. ( 1471) ). the piecewise expression of /oA( r j|) follows from 
Eq. (HHD as 



CW) 



In 



arctan 



/r?(l-r ) 



In 



+ 7ri 



(rjj < 0) 
(0 < rj < 1) 
(1 < r 2 ) 



(51) 



Notice that this result does not depend on B r , reflecting the fact that the threshold condition 
g 2 = 4m 2 is independent of B r . 

In principle, one could compute the vacuum polarization tensor by an alternative method. 
By using fermion propagators which are decomposed into the Landau levels 
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2a|, the one- 



loop diagram can be calculated order-by-order with respect to the Landau levels. The out- 
come would however look rather complicated compared with our calculation in the proper- 
time method. Nevertheless, as far as the lowest Landau level is concerned, a simple calcu- 
lation leads to an analytic expression of the vacuum polarization tensor at this level. This 
calculation was recently carried out in Ref. 30(, and the result is exactly the same as above. 
Therefore, our analytic results correctly reproduce the lowest Landau level contribution, 
suggesting the legitimacy of our interpretation of the Landau levels. 

Note that the vacuum birefringence is specified by the external photon momentum and 
the magnitude of the external magnetic field. Then, we notice that the lowest Landau level 
contribution is not necessarily enough for the calculation of the dielectric constants and 
refractive indices of energetic photons. Figure H] shows variation of the threshold structure 
with increasing magnetic field B r . The vertical line corresponds to the longitudinal photon 
momentum squared r 2 = gjj/(4m 2 ), and curves show locations of the thresholds r 2 = sfp. 
We have shown only up to the 14th threshold. The high momentum region above the 14th 
threshold line (i.e., the shaded region) is filled with narrowly spaced many curves. This 
figure visualizes that there is certainly a region (i.e., low momentum and strong field) where 
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FIG. 4. Threshold structure of photon momentum as a function of B r . Curves are positions of 
thresholds rjj = s+. We have shown threshold lines from the lowest (horizontal line at rjj = 1) up 
to the 14th. With increasing B T , the lowest threshold is isolated from higher levels more distantly, 
and thus the lowest Landau level approximation is justified in a wider kinematical region below 
the second threshold. 

the lowest Landau level approximation is appropriate. However, if a photon momentum rjj 
approaches higher levels, the levels close to the photon momentum would provide dominant 
contributions and thus higher-level contributions have to be incorporated in calculations 
even in strong magnetic field limits. With the results of a series representation in Eq. (1391) . 
important contributions from individual Landau levels can be systematically taken into 
account in calculating the vacuum polarization tensor. 

C. Behavior of the coefficients \i a * the thresholds 

Let us carefully see the behavior of the function at the boundaries = s±, which 

should be essentially the same as those of the coefficients Xi- First of all, given the explicit 
form of I^{rT) in Eq. (JUJ), one might expect that singularities could appear at both of the 
boundaries due to the inverse-square- root. But in fact, one finds that I^ir'i) is smooth 
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at the lower boundary rjj = s e ™ while it has a discontinuity with divergence at the higher 
boundary rjj = s+ (see Appendix [Hi for more details). 

As rjj approaches the lower boundary rjj = s!r from below and above, one finds that the 



values of i^(r 2 ) are finite and coincide with each other: 



lim i? A (rf) = lim J&(rJ} 



r?->sl™-0 rf.^-s^+O C — CL 



2p(£,n;B T ), (52) 



where a limiting value is given by p(i,n;B r ) = {(1 + 2£B r )(l + 2(£ + n)5 r )} -1 / 2 . Note 
that this is a real number as we commented before. On the other hand, if rjj goes beyond 
the higher boundary rjj = s+ (threshold for the photon decay), the function /^(rjj) has 
an imaginary part. Then we find that both the real and imaginary parts of i^(rjj) show 
singular behavior at the higher boundary rjj = sip. The real part goes 



A(r 2 



lim Re [/; A (rj 



2' 

) = UII1 ^=^= 

r},— >s^ — ^ 

9 



lim 7" A (rjf) = lim - " = +oo , (53) 



r^s^+O " c — a 



2p(£,n;B I ), (54) 



where the coefficient A (rjj) = 7r(rjj — s^ 1 ) -1 / 2 takes a finite positive value at the threshold, 
A(slp) = up? (£,n; B r ). On the other hand, the imaginary part diverges as rjj approaches 
r 2 = s+ from above, while it is just zero from below: 

lim Im^r 2 )] =0, (55) 

rf->s+™-0 

A(r 2 ) 

lim Im [L n A (r 2 )l = lim - " = +oo . (56) 



II ° + 

Such singular behavior appears only at the higher boundary rjj = sip for each number of i 
and n. The finite value p(£,n;B T ) and the coefficient A(rjj), as well as the position of the 
singularity, depend on the indices £, n, and the magnitude of the magnetic field B r , but the 
global structure of the function i^(rjj) is common for any £, n and B T . This is shown in 
Fig. El 

Recall that the coefficients x% are functions of -^.(r 2 ) and thus they will show the same 
threshold behavior. Since Xi are obtained after the summation over all the Landau levels 
specified by i and n + £, and each /^(rjj) has a singular point for each n and £, Xi have 
divergences at infinitely many thresholds. In Fig. |6l we show the real and imaginary parts 
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FIG. 5. The function is plotted against r| for different values of n and B r = 1. Blue (red) 

solid and dashed lines are real (imaginary) parts of ^^\( r y ) for n = and 1, respectively. Filled 
squares (circles) on the lines are lower (higher) thresholds r| = sf? (r| = s+ ). 

of Xi as functions of rjj for fi r = 10. Among infinitely many terms, we have summed 
from the first term up to the 14th term. Spikes correspond to the thresholds for different 
values of £ and n. For example, the first and second thresholds appear at rjj = s*} = 1 
and rjj = = (1 + y/l + 2i? r ) 2 /4 ~ 8 for B T = 10. It is important to notice that 
these divergences are harmless in the dielectric constants e±, ey and the refractive indices 
n±, n|| . This is easily understood from the explicit representation of the dielectric constants, 
Eqs. ( )23|) and (I24I) . First, the function enters all the coefficients Xi — 0, 1,2), and 

thus three of the coefficients have divergences of the same order at exactly the same momenta. 
Second, these coefficients appear both in the denominators and numerators of Eqs. ( 123|) and 
( )24|) . Therefore, the singularities of /^(rjf) bring the same divergences in the denominators 
and numerators, and are canceled to give finite values of the dielectric constants and the 
refractive indices. 

We also notice that magnitudes of the dielectric constants, in the vicinity of a threshold 
specified by £ and n, are governed by the prefactors of the terms proportional to /^(rjj) 
in because they dominate all the other finite terms including the terms proportional to 
/p' A (rj|) with indices n' different from £, n. 
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FIG. 6. Real and imaginary parts of the coefficient xi as functions of r?, for B T = 10. Shown 
are the results obtained with the first 14th terms. While thin colored lines indicate individual 
contributions from Landau levels shown with the same colors in Fig. HI a thick line indicates a 
resultant quantity after the summation. 
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V. SUMMARY AND PROSPECTS 



In this paper, we have presented analytic representation of the vacuum polarization tensor 
in a strong magnetic field, which is an indispensable quantity for the description of photon's 
birefringence. This was made possible by the use of infinite series expansion applied to 
the exponentiated trigonometric functions, whose integer indices have been later physically 
interpreted as the Landau levels of a fermion-antifermion pair. The infinite sum corresponds 
to the contributions from all the Landau levels to the vacuum polarization tensor. We 
have carefully investigated the singularity structure of the vacuum polarization tensor at an 
infinite number of the thresholds beyond which a photon decays into a fermion-antifermion 
pair. 

Having obtained the analytic expressions of the vacuum polarization tensor, we can pro- 
ceed to computation of the dielectric constants ejn and refractive indices n_ui by using 
Eqs. (123)) and (I24I) . However, this does not go straightforwardly as we have already empha- 
sized in the last paragraph of Sec. IIII1 We need to solve Eqs. (1231 and ( 1241) self-consistently 



with respect to the dielectric constants. In Ref. 18[, such self-consistent solutions were nu- 
merically obtained when the photon momentum is smaller than the lowest threshold rj? < 1. 
The result shows that one of the refractive indices starts to deviate from unity as the mag- 
netic field strength increases beyond the critical one B r > 1, and approaches a limiting value, 
depending on propagation angle, in the extremely str ong field limit B T 3> 1. Our calcula- 



tion of course precisely reproduces the results of Ref. 18J, and can go beyond it. However 



we have to be very careful when we investigate a region where the deviation of the dielec 



trie constant becomes sizeable. In the next paper [22J, we will inspect the self-consistent 
solutions when the photon momentum rjj increases across the threshold. 

The general expression obtained in this work beyond the threshold has opened new pos- 
sibilities to study a wide variety of phenomena which contain quite different scales. Recall 
the threshold structure shown in Fig. HI and note that the kinematical region discussed in 
covers a small region below the first threshold (horizontal line at rj, = 1). How- 



Ref. 



ever, if one wants to apply to relativistic heavy-ion collision, photons generally have large 
momenta, and thus we have to take into account the effects of many thresholds. Also, 
in application to the physics in magnetars, one would be interested in the interplay and 



competition among photon splitting, vacuum birefringence and photon decay 
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32|. To 



investigate these intriguing phenomena, we need to accurately describe the structure of the 



221 ]. we will 



photon spectrum with the thresholds taken into account. In the next paper 
discuss how and to what extent the dielectric constants and the refractive indices are in- 
fluenced by the singular behavior of the vacuum polarization tensor at the threshold. We 
will show that the location of the threshold itself has to be obtained in a self-consistent 
way. In particular, we will investigate behavior of refractive indices incorporating effects of 
thresholds, which can be studied within the lowest Landau level approximation. 
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Appendix A: Fermion propagator in a constant magnetic field 

We describe here how to obtain the dressed fermion propagator in a constant magnetic 
field shown in Fig. [TJ As mentioned in the text, we employ Schwinger's proper-time method 



20 
t: 
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24J, which yields the integral representation (JSJ) with respect to the "proper-time" 



1 f°° r i 

G(p\A d ) =i(f- e4d + m) x - / dr j*{<t-*4*) a -<r» a -*)} ? 

i> Jo 

where is the Fourier transformation of A^(x) and a small imaginary value — ie, originated 
from the prescription of the boundary condition, makes the integral convergent. An advan- 
tage of this expression is that the terms containing the gauge field appear in positive powers, 
while in the original definition the gauge field appears in the denominator. Below, we 
analyse the following integral part: 

i r°° 

A{p\A cl ) = - df e^W-^-i^-ie)} ) ( A1 ) 
* Jo 

which, equivalently to Eq. satisfies the following equation: 

{(jf - e4 cl ) 2 - m 2 } A(p|A cl ) = 1 . (A2) 

Note that the operator contains Dirac matrices, and thus both the equation and A(p|A c i) 
are understood as matrices with respect to spinor indices. 
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It is convenient to perform calculation in the Fock-Schwinger gauge for the external field, 
x^A^x) = 0, which allows one to represent the gauge field in terms of the field-strength 
tensor as = — \F ilv x v when the field strength is constant. An elementary calculation 

yields the expression for Eq. ( 1A2I) in terms of the field strength: 



p 2 - k 2 + — — — F^F, 



e d „,„,„ d 



4 dpv dp a 
where we have introduced k defined by 



A(p|Ai) = l , (A3) 



K 2 = m 2 - i € + -F^o^ , 

with the antisymmetric tensor, = f [7 M , 7 I/ ]- 

By solving Eq. ( 1A3I) . one can derive another expression alternative to Eq. ( 1A1I) . which is 
favorably written with respect to the field-strength tensor. Namely, we find 

A(p\A cl ) = - df e -^ 2 T+iP^(+)p„+Y(T) 
* Jo 



where we have introduced two functions X^ u {r) and Y(r) defined by 

i(ei 

7b 



X^(f) = [(eFy 1 tanh(eFf)]^ = f rjf + t&n ^) ^ j (A5) 



Y(t) = -^Tr [ln{cosh(eFf)}] = - hi |cos(e5f)| . (A6) 

In this expression, we have suppressed contracted tensor indices, and written inverse field- 
strength tensor as [F~ l ) tiu . As shown above, two functions X^ u {f ) and Y(f) are especially 
simplified when there is only a magnetic field: The field-strength tensor is then given by 
F 21 = —F 12 = B with all the other elements vanishing. We find that X^ v shows different 
behavior in the longitudinal (||) and transverse (_L) planes with respect to the external 
field, which are, respectively, specified by the metric tensors r]^ v = diag(l, 0, 0, — 1) and 
7?f = diag(0, -1,-1,0). 

Substituting the above expressions ( 1A4|) - ( 1A6I) into Eqs. (JSJ), one obtains the following 
integral representation of the dressed fermion propagator: 

POO 

G(p\A cl ) = df f- e^F^X ua (f)p a + m e -i^+ip^{r) Pv +Y(r) _ ^ 
Jo L J 

Recall that n 2 contains the gamma matrices, and thus we have to take into account the 

order among them when we carry out the trace of spinor indices. Although integration 

over the proper-time f has not been performed, this compact form allows us to efficiently 
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simplify awkward expressions in applying to resummed diagram calculations which include 
the vacuum polarization tensor shown in Fig. [2j Owing to the Fock-Schwinger gauge, a 
configuration of the external field is manifestly specified in terms of the field-strength. 



Appendix B: Renormalization of the vacuum polarization tensor 

Below we first discuss renormalization of the polarization tensor in the ordinary vacuum, 
and then consider the case with the magnetic field. 

Our framework for nonzero magnetic field background, of course, contains the case of 
the ordinary vacuum as a limit of vanishing magnetic field. This is explicitly checked in 
the expressions fl9]) - (TTTT) which have not been expanded with respect to the Landau levels. 
Indeed, in a vanishing magnetic field limit B T = 0, the scalar functions Xi = 0, 1, 2) in 
Eqs. (J9]) - (jTTj) reduce to the results in the ordinary vacuum. This is easily seen as follows. 
Since r x and T 2 vanish in the vanishing field limit B T — > with r and P being fixed! 14 ! the 
scalar functions Xi an d X2 also vanish in this limit: 

Xi = X2 = . 



On the other hand, Xo reduces to the following form 



xT(r 2 ) = —J d/3 J dr —JLe-^ , (Bl) 

where we have defined a phase <f)o(r 2 ) = 1 — (1 — /3 2 )r 2 with a normalized squared momentum 
r 2 = q 2 /(4m 2 ). We have introduced a cut-off 5 (0 < 5 1) at the lower boundary of the 
integral with respect to r, because, without regularization 5 — > 0, there arises a logarithmic 
divergence in this integral. In terms of Xo aC ; ^ e vacuum polarization tensor in the ordinary 
vacuum Ilo^g 2 ), obtained as the vanishing field limit, is given by 

nn?) = -xo a V)iT • (B2) 

When the squared momentum r 2 is less than the threshold^!. < r 2 < 1, the phase 4>q 
becomes positive definite, and then the integral with respect to r in Eq. (1B1I) is calculated 

as 



xS ac = ^ / dp(i-p 2 )r(o,<i> 6) 



14 



One should take this limit after pulling back the proper-time r to the original dimcnsionful form, r 



|eB|r, because the dimensionless form varies as the field-strength varies. 
15 Expression in the other range is obtained by an analytic continuation afterward, if necessary 
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= ~ J d^{l - (5 2 ) { 7E + In 0o + In 5 + O^ 1 )} (B3) 

where, on the first line, we have changed the contour so that we can use the 'incomplete 
gamma function', T(s, t) = f°° x s ~ 1 e~ x dx, and, on the second line, we have expanded it with 
respect to the small cut-off scale S. The constant 7e denotes the Euler number. Note that 
a divergence in the limit 5 — > is isolated from the other finite terms. 

Adopting the conventional on-shell renormalization condition, we define a finite vacuum 
polarization tensor U.Q U (q) in the ordinary vacuum as, 

nn g ) = nn?) - n(T(o) • (B4) 

From the above on-shell condition and the expression (1B3|) for Xo ac ? ^ i s natural to introduce 
a finite scalar function Xo a °( r2 ) by 

Xir(r 2 ) = xS*(r 2 ) - C(0) (B5) 
= ^/W-^) to { 4m2 _^ w }. (B6) 

Namely, the finite vacuum polarization tensor is given by Xo a °( r2 ) : 

flT(ff) = -*r(r 2 )^r, (B7) 



which agrees with the conventional result (see, e.g, Eq. (7.91) in Ref. |23|). An integral in 
Eq. ( 1B6I) with respect to (3 is carried out analytically, and then we find an analytic form, 
which allows for an analytic continuation to the regimes r 2 < and 1 < r 2 , as 

a 



D7T 



3^\^r 2 /]l^ 2 / w 2 



(B8) 



Now let us turn to the case with magnetic fields. As mentioned below Eqs. (l9~))-(llip. the 
vacuum polarization tensor in the presence of the external magnetic field contains exactly 
the same logarithmic divergence as the one shown in Eq. flB3p . Therefore, we remove the 
logarithmic divergence by imposing a similar on-shell renormalization condition, 

U^(q h q ±] B t ) = rC( % qr, BJ - Tl^(q 2 = 0) , (B9) 

so that the finite vacuum polarization tensor TL%£(q\\, qj_; B T ) vanishes at simultaneous limits 
of q 2 = 0, q\ = and B r = 0. In terms of the scalar functions \ii the renormalization 
condition is written as 

n^(g,|,^; Br) = - {xo(rf, r\- B r ) - x ™(0)} ^ " E ^ ■ ( B1 °) 



i=l,2 
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To proceed further, we have to regularize the divergent scalar functions xo and Xo ac - 

We show regular izat ion methods which are involved in our numerical 22| and analytic 
calculations. First, in case of numerical calculation, we regularize xo an d Xo ac by inserting 
common cut-off scales into these functions. Then, the scalar function xo hi the presence of 
the external magnetic field is regularized in the same manner as in Eq. (IBip . 



Xo(rLr 2 ± ;B T 



a f 1 



sin r 



(BH) 



rf) + 



where we have put all the phase appearing in Eq. ([9]) together as 0(r|,r^) = 
(wcos(/3r) — i] cos t)/t. However, unlike the case in the ordinary vacuum (lB3p . a complex 
structure of xo prevents us from isolating the divergence by an analytic calculation. We 
therefore numerically evaluate a convergent integral, 



a 
Air 



1 PCX} 

d(3 \ dr 
i Js 



r Hr) 

sinhr 



/3 2 



-0o(O)r 



T 



(B12) 



where we have rotated the integral contour downward in the complex r-plane assuming a 
momentum regime rj? < (1 + y/1 + 25 r ) 2 /4 so that the integrand vanishes on the large arc. 



We then obtain a finite part of the scalar function as, 

Urlrl-,B T )= XQ {rlrl-B T )-xT^) 
= J (rlrl;B r )+K (6) 



(B13) 
(B14) 



We have put a remainder as 



a 
in 



1 i-S 

d/3 / dr 
i Jo 



r (-ir) 



r 2 ± )i 



I"/? 2 



>o(0)r 



(B15) 



sinh r 

which gives a cut-off dependence on the right-hand side of Eq. flB14l) . Note that, if we could 
remove the cut-off scale by taking a vanishing limit 5 — > at the end of calculation, the 
renormalized scalar function Xo( r y ? r ±] B r ) should not depend on the cut-off scale as long 
as we insert common cut-off scales into Xo an d Xo a °- However, we cannot take this limit in 
numerical calculations, and thus it should be estimated how much the renormalized scalar 
function Xo depends on the cut-off scale 5. For a small r, the integrand is composed of the 
following functions, 

Tui-tr) 1-/3 2 , 1 



~ 1 — 

sinh r r 6 

)(^,ri)~0 o (r 2 ) + 0(r 3 ) 



1 - f3 2 ) 2 T + 0(r d 



(B16) 
(B17) 
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and thus the remainder is evaluated as, 



a 

AlT 

4a 



dp / dr [-(1 - p 2 ) {Mr 2 ) ~ <M0)} + <D(t)] 



15tt 



r 2 5 + (D{5 2 



(B18) 
(B19) 



The above estimate shows that the remainder TZq(5) is linearly suppressed as we take a small 
cut-off scale 6 in numerical calculation. 

There is an improved method to suppress the cut-off scale dependence, which we have 
adopted in our numerical calculation. Note that the term independent of r in Eq. (1B18j) 
vanishes, if two phases <fto are evaluated at the same momentum. We rearrange Eq. (IB 13j) 
to eliminate this constant term by using Eq. (1B5P as 

Urlr 2 ± ;B r ) = { Xo (rJ, r\- B r ) - xT(r 2 )} + ^(r 2 ) . (B20) 

The renormalized scalar function is then obtained as, 



Xo(rJ, ri; B t ) = J(rJ , r\; B t ) + xT(r 2 ) + K{8) , 
with a convergent integral J(rn,r^; B r ) and a remainder TZ(S) given by 



a 



J(rtri;B T ) = — / dp 



K{5) 



47T 

a 
4ir 



dr 



dp / dr 



To(-ir) 
sinh r 

r (-ir) 



(rlrDr 



1-/3 2 



-<)4o(r 2 )T 



(B21) 

(B22) 
(B23) 



sinh t r 
We already have an analytic expression for the renormalized scalar function Xo a °( r2 ) m the 
ordinary vacuum. In this method, the remainder TZ(5) is evaluated as, 



a 

An 
a 



dp / dr 



(l-^(H)--(l-^Y \t + 0(t 2 ) 



6 



1807T 



(11-12^)^ + 0(5 



(B24) 
(B25) 



By using this method, the cut-off dependence is improved from a linear dependence Hq{8) ~ 
0{8 X ) to quadratic one K{$) ~ C(5 2 ). 

Next, we proceed to the case of our analytic calculation. In the expansion method 
resulting in the double infinite sum with respect to the Landau levels, we do not have to 
insert a cut-off scale to perform term-by-term integrals, because these contributions, from the 



10 



An equivalent expression to Eq. (|B20p has been shown in Rcf . 19( , and might be adopted to numerical 



computation in Rcf. 



1 



18[. However, it has not been clear how it improves the naive subtraction (|B13() . 
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discretized levels in the transverse momentum, does not diverge individually. However, by 
taking all-order summation, we show that the infinite sum leads to a logarithmic divergence 
which should be identified with the logarithmic divergence we have seen above. 

We examine the double infinite series at vanishing momenta = and r\ = 0, because 
the ultraviolet divergence itself does not depend on external momentum scales. At these 
momenta, the series coefficient C" and its multiplicative form in xo have simple expressions, 

, 1 (n = 0) 

C £ n (0) = { y ' (B26) 
(n > 1) 



n 



1+1 (n = l) 

-Cn0)=< (B27) 
V |^ (n > 2) , 

where we have used L®(0) = 1 and Lj(0) = I + 1. Note that the expression of xo should be 
understood that any term being proportional to x C™(0) at n = is absent, because 
such term never exists in the expansion flE3j) . Then, summation of the coefficients in \2 (see 
Eqs. ( )43|) and ( 144)) ) are also easily performed as 

i-\ 

^ 0(1) (0) = -8^(£-A){-Q°(0)}=4^ + l) , (B28) 

A=0 
t-l 

D; (1) (0) = -8^-A)C°(0) = -U{1+1) , (B29) 

A=0 

£-2 

Dj (1) (0) = -8J3^-A-1){-C?(0)}=^-1) . (B30) 

A=0 

Therefore, the infinite series with respect to n terminates up to the first two terms. At 
vanishing momenta and for n — 0, 1, the functions i^ n (0), (j"(0) and H™(0) have the following 
forms: 

= TvIts- , (B31) 



F}(0) = ~E\, (B32) 
G K0) = ~ [{1 + (2f + l)B r } E] + 2B T ] , (B33) 



if °(0) = ^— . (B34) 

By using these explicit expressions, we show that the scalar coefficient xo an d differences 
Xi~ Xo an d X2~ Xo contain the same logarithmic divergence. Introducing Xi = Xi~ (1 — ^o)Xo ; 
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we have 



X t (rf l =rl = 0;B r ) 



An 



u 



(B35) 



where we separated an infinite series with respect to £ so that it begins from £ = 1. Explicit 
expressions of terms are given by 

f 2{F^-Gl) (i = 0) 

F ° + F° - + #°) (< = 1) 
(z = 2) 



Ui(B T ) 



(B36) 



and 



(B37) 



2{F}-(2£ + l)G}} (i = 0) 

F° + - (flj + £T° +1 ) (i = l) . 
^ + l){(F°-2F/) + F° +1 } (i = 2) 

Here, we suppressed the vanishing arguments of F™(0), G™(0) and Hg(0). To see a behavior 
of the summation taken up to a large £, we insert a cut-off at £ = L assuming an infinite 
limit L — > oo. The infinite sum is then evaluated by being replaced by an integral as, 



Xi(rf\=rl = 0;B T 



^\ Ui (B r )+ hm LV^(l;B r )) 

All L^oo L 

Ui(B T )+ lim L [ dxWi(Lx;B T )\ . (B38) 



aB 
An 



The above integrals are carried out analytically, and then expanded for a large value of L 



as, 



Xi{r\\ 



0; B t ) = \n(2B t L) + V c k {B T )U 



(B39) 



fc=0 



where series coefficients are denoted by Ck{B T ). The above expansion series does not have 
any positive power of L, and thus the leading term for large L (with finite B T ) is given by a 
divergent logarithm. Note also that a coefficient of the divergent term is independent of B r , 
implying the same divergence as in the ordinary vacuum. While the scalar functions Xi an d 
X2 do not contain any divergence owing to a cancellation between the logarithmic terms in 
the functions Xii the other one xo contains the logarithmic divergence found in Eq. (IB39I) . 
These results are consistent with the observation below Eq. OU) - fITTjl . Therefore, we might 
be allowed to identify this divergent logarithm, appearing from the infinite sum, with the one 
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FIG. 7. A modified photon propagator 

in Eq. ( 1B3I) appearing in the limit 8 — > 0. Comparing them up to the divergent logarithmic 
term as 

f InOr 1 ) = £ ln(2B r L) , (B40) 

37T 37T 

we then find a relation between the cut-off scales, 

cT 1 = 25 r L . (B41) 

According to this relation (lB4ip . the ultraviolet logarithmic divergence in the presence of 
a magnetic field is subtracted by the same logarithmic divergence in the ordinary vacuum. 
Namely, the divergent contribution in a sum up to £ = L ^> 1 might require to be subtracted 
by an integral (1B1I) evaluated with a cut-off scale 5 = (2Z? r L) _1 , and a finite part in our 
prescription is defined as a remainder of this subtraction. 



Appendix C: Dispersion relations from modified photon propagator 

Here we demonstrate that one can obtain Eq. (|2ip (that immediately yields the dispersion 
relations ( 123]) and (124|) ) from pole positions of a modified photon propagator after taking 
into account the interactions with the external magnetic field. 

First of all, it should be noticed that insersion of the self-energy 11^ into a photon 
propagator corresponds to resumming the dressed-fermion one-loop diagrams as shown in 
Fig. 7. There, the dressed-fermion one-loops are linked by the free photon propagators 
D^ v {q 2 ) given by 



(Cl) 



A full photon propagator in this approximation can be expressed as 

oo 

LH ? )= J>n?), (C2) 



n=0 
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where the n-th contribution corresponds to the diagram having n rings, and n = is the 
free propagator Dq V . To compute the n-th order contribution D£ u (q 2 ), we use the following 
identities among the projection operators P? v : 

PZ P P ip v = P? P P* p v = q 2 Pr (* = 0,1,2), 

p^p,; = qfpr , p% p p 2p v = <&? , *r *y = p% p p Xp v = o . (C3) 

Since, as shown by the above relations, the products of two projection operators always 
result in zero or one of the projection operators, we are able to expand the n-th order 
diagram into the same form as the self-energy with respect to the projection operators: 

= C (n)p r + C (n) pr + c (n )pr (Q4) 

For n > 1, we can drop those terms proportional to q^q u in the free propagator fICll) . 
because they vanish when contracted with any of the projection operators in the self-energy 
as qfiPf = (i = 0, 1, 2). The coefficient scalar functions are then calculated as 

Ci n) = (-)x(-^) (?T-\ (C5) 

n—r 

c!"' = f ^ ) x T.nC r I -§) ( -a ) (fy 




g 2 7 

(XoA , (C6) 

C^ = Ci (n) (1^2, , (C7) 

where the last line is meant for replacements of subscripts in ([ n ^ as indicated. To carry out 

the summations in £i and Q , we used a formula for the binomial expansion. Substituting 
Eq. fjC4j) with the coefficients (1C5j) - (1C7I) into Eq. (1C2|) . we obtain the resummed photon 
propagator as 



iD^(q) 



i PT , ^i(g 2 ^y) , Mg 2 ,gl) iff" , & ?<t 

'X + xo)<? 2 <? 2 (i + xo)q 2 + xiqf qf (i + xo)<? 2 + X2<?i ?i <? 2 <? 2 

^0)^0) , n w n w , ^2)^2) ^ u „ rrsl 

1 + Xo )g 2 (l + xo)<? 2 + x 2 <?i (l + xo)q 2 + xig,f <? 2 (3) (3) ' 1 J 
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We have introduced a|j(g 2 ,g 2 ) = -gjxi/{g 2 (l + Xo)} and a ± (q 2 ,q 2 ± ) = -q 2 ± X2/{q 2 (l + Xo)}, 
and utilized relations shown in Eq. f|T7|) to get a 'diagonalized' form on the second line. 
Contracting the propagator (JUS]) with eigenvectors as, n^D^n^ {i = 0,1,2,3), we find 
that the pole position of each mode exactly agrees with what has been obtained in Eq. ( 12"T]) 
from diagonalization of the modified Maxwell equation. 



Appendix D: Physical modes of a photon in a magnetic field 

In this Appendix, we show that, among four polarization modes in Eq. (12"T|) . only two 
modes, having nontrivial dispersions, are physical. We consider electric and magnetic fields 
of the dynamical photon field (120]) whose polarizations are specified by the eigenvectors 
17 (\) (-^ = 0, • • • ,3) defined in Eq. (TT9T) . Explicitly, the Fourier modes of the electric field 
£ l = —F 0% and magnetic field B % = \e l ^ k F^ of a propagating photon in Eq. (|20]) are, 
respectively, given by 



& {x) (q) = iN (e^'V^A)) e^<^ . 



(Dl) 
(D2) 



Let us first consider the polarization modes A = and 3 (see Eq. ( 1T6]) ). both of which give 
the ordinary dispersion relation of the massless type u 2 = \q 2 \. Substituting the definitions 
of 71"^ and 7r^ into the above, one immediately finds 



'(3) 



£(0) oc 



Qy 



B 



(0) 



for A = 



(D3) 



and 



5 ( 3) =0, 4) = 0, for A = 3, (D4) 

where the dispersion relation q 2 = has been used for the A = mode. Namely, the two 
modes A = 0, 3 do not induce excitation of the electromagnetic fields, which implies that 
they are unphysical. In fact, these modes are related to redundant gauge degrees of freedom. 

In contrast, the remaining two modes A = 1,2 provide nonvanishing electromagnetic 
fields. First of all, one can easily verify (without using the dispersion relations) that the 
electric field and the magnetic field are orthogonal to each other for each mode: S^\) • £>(A) = 
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(A = 1,2). Next, let us employ the same coordinates as those used in deriving Eqs. (|23|) 
and (I24p . Then the electromagnetic field of the A = 1 mode (with momentum g M ) reads 



x) = tuN 



v / 



where the common factor ^±(t, x) is defined by 

ty±(t, x) = exp{ - iu(t - n v l al q ■ x)} exp{ 



' cos 9^ 


\ sin 9 j 



imae * ~i 

u>n ± q ■ Xj . 



(D5) 



(D6) 



Note that the factor includes real and imaginary parts of the refractive index n±, which 
respectively represent oscillation and decay behavior of the electromagnetic field. Thus, 
they give a phase velocity and damping coefficient, respectively. A unit vector q denotes 
a three-dimensional vector directed in the photon momentum q. As in a photon in the 
ordinary vacuum, these fields are transverse to the photon momentum: q ■ = q ■ B(i) = 0. 
The electric field induced in this mode is oriented to the transverse direction with respect 
to the external magnetic field, and thus does not have parallel component. 

Similarly, the electromagnetic field of the last remaining mode A = 2 is given by 

( e,| (f ) cos6» \ / ^ 

V\\(t,x), B {2) =coN 



C2) 



(0) 



V 



sin 9 



7T 

'2 



/ 



*\\(t,x).(D7) 



V / 



where e»(9) denotes the dielectric constant when the photon momentum is directed to the 
zenith angle 9. While the magnetic field is transverse to the photon momentum q ■ B^) — 
0, we notice that the electric field in general has an overlap with the photon momentum 
as, q- Bp) 



uN 



gr {ey (l) — l} sin 9 cos 9, depending on the propagation angle 9 with 



Hl(f) 



respect to the external magnetic field. 

The space-time dependence of the field is again represented by a common factor ^f\\(t, x), 
and it describes oscillation and decay behavior of the field. However, it contains the refractive 
index n» (instead of n±). Namely, the factor is explicitly written as 



\&ll(t, x) = exp { - iu(t - n|j cal q ■ x)} exp { - 



a;nj| mag q ■ x } . 



(D8) 



Recall that the refractive index nn is in general different from n± in *$?±(t, x). Therefore, the 
two physical modes A = 1,2 have distinct phase velocities and decay rates (into a fermion- 
antifermion pair). The electric field shown in Eq. (1D7|) contains a component oscillating 
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along the external field in which direction a charged particle has continuous spectrum. This 
is in contrast to the electric field of the other physical mode. Since the external magnetic 
field provides anisotropic spectra of fermions in the Dirac sea, oscillating electric fields, 
and £{2)1 of an incident photon induce anisotropic polarizations. 



Appendix E: Partial wave decomposition 

We expand exponentiated trigonometric functions in Eqs. (Q - (fTTj) by using the Jacobi- 



Anger expansion that is well-known as the partial-wave decomposition 26]: 



exp {— iu cos(/3r)} = ^^(2 — 5 n o)(—i) n J n (u) cos(n(3r) 

n=0 

oo 

= ^(2 - S n0 )I n (-iu) cos(n(3r) 



n=0 

oo 



J> - 6 n0 )I n (-iu) e m ^, (El) 



n=0 

d 

cos(/3r) exp {— iu cos(/3r)} = i— exp {— iu cos(/3r)} 



du 

-(2 - 5 n0 ) {I n+1 (-iu) + cos(n/3r) 



2 

n=0 



oo 



^Y.V 2 - <W { J «+i(-^) + e m/3T , (E2) 



2 

n=0 
1 9 

sin(/3r) exp {— iu cos(/3r)} = — -— — exp {— iu cos(/3r)} 

oo 

2 



^ n(2 - S n0 )I n (-iu) sm(n/3r) 

n=0 

-, oo 

-Vn(2-y4(-m)e*, (E3) 



n=0 

oo 

1 ■ 



M n=0 



where J n (.2) and are the Bessel and modified Bessel functions of the first kind related 
as, I n (z) = (—i) n J n (iz) . Noting that the integrands in Eq. (|2"5|) are even function of /3, we 
have replaced the trigonometric functions by exponential factors on the right-hand sides of 
the arrows, i.e., cos(n/3r) — > e ml3r and sm(n(3r) — > — ze m/3r . 
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Appendix F: Limiting form of the coefficient Cf{rj) 

An approximate form of C™(rj) for large values of £ 3> 1 is found by the use of a limiting 
form of the associated Laguerre polynomial and Stirling's formula for the factorial £\. When 
£ is large enough with fixed n and 77, a limit ing form of the associated Laguerre polynomial 



L™(r]) is expressed as (see (13.5.14) in Ref. [331 ]) 



Lftifl ~ -1= ( " V * ei cos (Fl) 



* ? fo) = 2 A M* + ^)-|(n + ^) . (F2) 



Also, we use Stirling's formula for large £ (see (6.1.38) in Ref. [33|). 

ei~ y /2^ e - i e i+ * . (F3) 

Then, the ratio of factorials in C"(r/) for large £ ^> 1 reads, 

£\ ( £ \ e+n+ ^ 



(£ + n)\ \£ + n^ 

^ e -^£- n , (F4) 



where we have used {£/(£ + n)} e+n+1/2 = e -(e+n+i/2)Hi+n/£) _ e - n {i+(n+i)/(2^)} to obtain the 
final expression. By using Eqs. (IFip . (1F2|) and flF4j) . we find a limiting form of C™{rf) as, 

C?fo) ~ -^T cos2 ( £ » X ) • ( F5 ) 

The amplitude of C™{rj) is bounded by a pref actor of the squared cosine which decays by an 
inverse-square-root as £ becomes large (see Fig. [3]) . 



Appendix G: Imaginary part of the function i]X( r f) 

In this Appendix, we carefully investigate analytic properties of the function If&{rV) = 
Ii + {ff) — I2-( r f) ( see Eq. ( 135]) ) which is the only one possible source of the imaginary 
contribution and thus determines the essential behavior of the vacuum polarization tensor. 
In the following, we consider the function I^{vT) with the indices £ and n being fixed. 

As mentioned in the text, the property of If^{rT) is specified by the discriminant D 
defined by Eq. (T4B1) . Using the solutions rj? = s± to the equation V = (see Eq. (14T|) ). we 
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consider three kinematical regions (for notational simplicity, we suppress the indices i and 
n in s±): 

(I) rj < s_ (V < 0) 

(II) s_ < r\ < s + (V > 0) (Gl) 

(III) s + <r\ (V < 0) . 

In region (II) where D > 0, -^7\(fjj) is obviously a real- valued function of r 2 , and thus 
there arises no imaginary part. In regions (I) and (III), the argument of the cotangent 
in If ± (r 2 ) becomes a pure imaginary number. The expression of I™±(r 2 ) in these regions 
is available via analytic continuation of the function in region (I). With an infinitesimal 
imaginary displacement r 2 + ie, analytic continuation into the regions (I) and (III) leads to 
the following form valid in both regions: 

— 2i 

Ie±( r f\ + = = arctan {— i(w± + ieK±)} 

\/b 2 — Aac 

= = {ln(l — w± — ien±) — ln(l + w± + ien±)} (G2) 

s/b 2 — Aac 

where we have introduced two real- valued functions, w±(r 2 ) — (b ± 2a)/\/b 2 — Aac and 
K ±( r jj) = b =F 2a. The latter functions k± are leading-order coefficients in an expansion of 
w±(r 2 + ie) with respect to the infinitesimal displacement. Owing to the expression in terms 
of the logarithm, we notice that the arctangent has an imaginary part when the absolute 
values of w± are larger than one, \w±\ > 1, and that its sign can be either positive or negative 
depending on the sign of k±. 

First, we show in what region I™±(r 2 ) has an imaginary part. Comparing the denominator 

2 

and numerator of w± as (b ± 2a) 2 — y/b 2 — Aac = Ar 2 {1 + (2£ + n =p n)B T }, we obtain 
inequalities \w±\ > 1 and \w±\ < 1, when the longitudinal momentum squared is positive 
r 2 > and negative r 2 < 0, respectively. Therefore, I™±(r 2 ) has an imaginary part only in 
the positive regime, r 2 > 0. 

Next, we should notice that the difference i]X(r 2 ) = I™ + (r 2 ) — I™_(r 2 ) has an imaginary 
part only when the functions I™ ± {r 2 ) have imaginary parts with opposite signs: otherwise, 
they cancel each other. The difference Ii&(r 2 ) retains an imaginary part when the arguments 
of the arctangent, w± + in±, are on the opposite sides of the branch cut. Thus, we shall 
examine the signs of w± and k± in more detail, which simply follow from those of b + 2a 
and b — 2a according to their definition below Eq. flG2j) . 
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In the following, the signs of the arguments are found in the regions, < rjj < s_ 
and s + < rjj (0 < s_ < s+), separately. First, we find an obvious inequality, b — 2a = 
—{nB x + 2rjj) < 0, as long as rjj > 0. Thus, the signs of sgn(w_) and sgn(ft + ) immediately 
follows in the both regions as 

sgn(u?„) < and sgn(re + ) < (0 < rjj) . (G3) 

On the other hand, the signs of the other functions depend on the momentum regime. To 
see this, we show two relevant inequalities valid in the each regime, 



b + 2a < -nB r + 2s_ < (1 + 2£B r ) - y/l + 2£B T =0 (0 < rjj < s_) , (G4) 
b + 2a > -nB t + 2s + > -nB T + - y/l + 2{£ + n)5 r 2 > (s+ < rjj) . (G5) 

Led by these inequalities, we obtain the signs of the remaining two functions depending on 
the regime: 

sgn(u> + ) < and sgn(K_) < (0 < rjj < s_) , (G6) 
sgn(w; + ) > and sgn(«_) > (s + < rjj) . (G7) 

Remind that the branch cut of the arctangent is located on the imaginary axis. From 
Eqs. (1G3|) and (1G6j) . we find that the imaginary parts of i]±(rjj) have the same signs in 



the regime < r^ < s_, and that the function i^(rjj) does not have any imaginary part 
in this regime. On the other hand, from Eqs. ( 1G3I) and (1G7I) . we find that the imaginary 
parts of It±{r^) have the opposite signs in the regime s + < rjj, and that the function i^(rjj) 
has an imaginary part in this regime. Therefore, we conclude that the difference /£X(rjj) 
retains an imaginary part only in the regime s + < rjj, namely, in region (III). A piecewise 
representation of /^(rjj) is shown in Eq. (j4"B|) . It has the same form for the real part in 
regions (I) and (HI), but has an imaginary part only in region (III). 

The imaginary part in the regime s+ < rjj is related to the analytic property of the 
vacuum polarization tensor in the complex r|-plane. It reflects that there are an infinite 
number of branch cuts initiating at the thresholds of photon decay, = sfp, specified by 
pairs of Landau-level indices I and £ + n, where a created fermion and antifermion pair has 
vanishing longitudinal momenta along the external magnetic field. Because the longitudinal 
momenta of a created fermion pair is continuous, each branch cut run continuously along 
the real axis into the asymptotic regime, ri? — > oo. 
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Appendix H: Limiting behavior of If^(r 2 ) at the boundaries 

In this Appendix, we discuss the limiting behavior of I^{r 2 ) m Eq. (148]) at the bound- 
aries, rj? = s± (we suppress the indices i and n for simplicity). We put a square root of the 



2 lim lUr\) 



discriminant as 5 = W \b 2 — Aac\ which vanishes at both of the boundaries. For infinitesi- 
mal 5, the logarithm and arctangent in I^{r 2 ) behave as log | | = + 0(5 3 ) and 
arctan{(6±2a)/<5} = ^■sgn(b±2a) — -^^5 + 0(5 3 ), respectively. Thus, the limiting behavior 
of /fA^jj) i s represented as 

' i + <W (1 -+ s ~ ~ 0) 

{sgn(6 + 2a) - sgn(6 - 2a)} TrcT 1 

(HI) 

+^ + 0{5 2 ) (rjj s_ + and s+ - 0) 

k 27TZO" 1 + ^ + 0(o 2 ) (rf S+ + 0) 
where we have used 6 2 — >■ 4ac (<5 — >■ 0) to obtain the zeroth-order term in the middle regime. 
Note that a negative-power term in the middle regime exists depending on a relative sign 
between the factors 6±2a. In Appendix [Gj we have already examined those signs just above 
Eq. (1G3P and in Eqs. (lG4p and (lG5p . They have the same signs as b + 2a < and b — 2a < 
in the vicinity of the lower boundary r 2 = s_ ± 0, while the opposite signs as b + 2a > 
and b — 2a < in the vicinity of the higher boundary r 2 = s + ± 0. [HI Therefore, -^(r 2 ) 
is finite in the limit r 2 — > s_ + 0, and, on the other hand, i2i( r f) ^ s divergent in the limit 
r 2 — >■ s + — 0. 

The limiting behaviors of I^{t 2 ) are summarized as follows. At r 2 = it is continuous, 
and thus has a limiting value I^{r 2 — > s_) = 2/(c — a)| r 2 =s _. On the other hand, /^\(^ 2 ) 
diverges to the positive infinity in the limit r 2 — >• s + — 0. The real part of I^ir 2 ) has 
a limiting value /^\(^ 2 — > s+ + 0) = 2/(c — a)| r , 2=s+ in the limit r 2 — )■ s + + 0, while the 
imaginary part diverges to the positive infinity in this limit. 
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